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Abstract. Using Fourier frequency resolved X-ray spectroscopy we study short term spectral variability in lumi- 
nous LMXBs. With RXFE/PCA observations of 4U1608-52 and GX340+0 on the horizontal/normal branch of 
the color-intensity diagram we show that aperiodic and quasiperiodic variability on ~ sec - msec time scales is 
caused primarily by variations of the luminosity of the boundary layer. The emission of the accretion disk is less 
variable on these time scales and its power density spectrum follows Pdisk(/) oc / law, contributing to observed 
flux variation at low frequencies and low energies only. The kHz QPOs have the same origin as variability at lower 
frequencies, i.e. independent of the nature of the "clock", the actual luminosity modulation takes place on the 
neutron star surface, 

The boundary layer spectrum remains nearly constant in the course of the luminosity variations and is represented 
to certain accuracy by the Fourier frequency resolved spectrum. In the considered range M ~ (0.1 — l)MBdd it 
depends weakly on the global mass accretion rate and in the limit M ~ M^dd is close to Wien spectrum with 
kT ~ 2.4 keV (in the distant observer's frame). The spectrum of the accretion disk emission is significantly softer 
and in the 3-20 keV range is reasonably well described by a relativistic disk model with a mass accretion rate 
consistent with the value inferred from the observed X-ray flux. 
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1. Introduction 

It is commonly accepted that in non-pulsating neutron 
star X-ray binaries the magnetic field of the neutron 
star is weak enough and the accretion disk can extend 
close to the surface of a neutron star. If the neutron 
star rotation frequency is smaller than the Kcplerian fre- 
quency at the inner edge of the disk, a boundary layer 
will be formed near the surface of the neutron star in 
which the acc reting matter decelerates and spreads over 
star's surface jSu nvaev fe Shakuralll986t lKluzniakLll 988: 
Ilnogamov fc Sunvaevl Il999t IPopham fc Sunvaevl l200ljh 
For a non-rotating neutron star, in Newtonian approxi- 
mation half of the energy release due to accretion would 
take place in the boundary /spread layer. The effects of 
the general relativity can increase this fraction, e.g. up to 
2/3 in the case of a neutron star with radius i?NS = 3i? £ 



(|Sunvaev fc Shakuralll98rilSibgatullin fc Sunvaevl |200C|) . 
Rotation of the neutron star and deviations of the space- 
time geometry from Schwarzschild metric further mod- 
ify the fraction of the energy released on the star's sur- 
face. Consequently, a luminous spectral component, cor- 
responding to the boundary layer emission is expected be 



Send offprint requests to: gilfanov@mpa-garching.mpg.de 



present in the X-ray spectrum of a neutron star X-ray 
binary 

X-ray observations of neutron star LMXBs in the 
high luminosity state reveal rather soft composite X- 
ray spectra. Based on theoretical expectations they are 
usually represented as a sum of two components at- 
tributed t o the optically thick emi ssion of the accre- 
tion disk ijShakura fc Sunvaevl Il973h and of the b ound- 
ary layer / neutron star surface ijMitsuda et all Il984l 
IWhite et allll98^1 . The spectra of these two components 
are rather similar to each other and decomposition of the 
X-ray emission into the boundary layer and accretion disk 
components is often ambiguous, especially when based 
on the spectral information alone. Not surprisingly, the 
best fit parameters derived from the data of different in- 
struments and, correspondingly, the inferred values of the 
physically meaningf ul quantities ar e often in contradiction 
to each other (e.g. iMitsuda et all Il984l iDi Salvo et all 
l200ll bone et all 12002^ . A robust and sufficiently model 
independent method of separating the boundary layer and 
disk emission is of interest. 

IMitsuda et all (|l984|) and IMitsuda fc Tanakal l)l986h 

studied pattern of spectral variability on the timescale of 
~ 10 3 sec in luminous LMXBs and found that the ob- 
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Fig. 1. The color-intensity diagram of GX340+0. Dashed 
line polygons show the regions at the Horizontal and upper 
half of the Normal Branch used for the frequency resolved 
analysis. 



served spectra can be represented as a sum of two compo- 
nents, having drastically different variability properties: a 
strongly variable ~ 2 keV nearly blackbody component 
and a stable softer component. They interpreted the hard 
and soft components as emission from the neutron star 
surface and from the optically thick accretion disk respec- 
tively and concluded that the optically thick disk is sta- 
ble on the timescales conside red. The latter conclus ion is 
in accord with the finding of IChurazov et all ((200lf> . who 
showed that the same is true in the high luminosity state 
of the black hole binary Cyg X-l. On time scales ^ 10 2 
sec the fractional amplitude of variations of the disk emis- 
sion is at least an order of magnitude lower than that of 
the hard Comptonized component. These results indicate 
that a low level of variability might be an intrinsic prop- 
erty of the optically thick accretion disk, independent of 
the nature of the compact object. 

As is well known (see Ivan der KlisL Il98fit 1200(1 for re- 
view) , aperiodic variability of X-ray flux from LMXBs can 
be broadly divided into two main phenomena - contin- 
uum noise and quasi-periodical oscillations (QPO) with 
frequencies r anging from several mHz to more than a thou- 
sand Hz (e.g.lHasinger fc van der KlisLll989Hyan der KlisL 
l200(ilRevnivtsev et allbOOlh . lMitsuda et alJl|l984 stud- 
ied the difference between the spectra averaged at differ- 
ent intensity levels - that restricted the range of accessi- 
ble time scales to <; 10 3 sec. In this paper we will exploit 
the technique of Fourier frequency resolved spectroscopy 
I Revnivtsev et all Il999l) to study spectral variability of 



Fig. 2. Power spectrum of GX340+0 at the upper part of 
the Horizontal Branch of the color-intensity diagram. 

luminous LMXBs on a broad range of time scales, includ- 
ing kHz QPO. 

As defined in lRevnivtsev et al] l|l999|) . the Fourier fre- 
quency resolved spectrum is the energy dependent rms 
amplitude in a selected frequency range, expressed in 
absolute (as opposi te to fractional) uni ts. A similar ap- 
proach was used bv lMendez et al to study the en- 
ergy spectrum of kHz oscillations in 4U0614+09. One of 
it's advantages over simple fractional rms-vs. -energy de- 
pendence is the possibility to use conventional (i.e. re- 
sponse folded) spectral approximations in order to de- 
scribe the energy dependence of aperiodic variability. One 
should keep in mind, however, that the interpretation of 
the frequency resolved spectra often is not straightfor- 
ward and might strongly depend on a priori assumptions. 
Nevertheless, several applications of this technique to vari- 
ability of black hole binaries gave mean ingful results (e.g. 
iRevnivtsev et all Il999t iGilfanov et afll200C|) . 

The simplest situation, when frequency resolved spec- 
tra can be easily interpreted is illustrated by the following 
example. Consider a two-component spectrum, in which 
one component is stable and the normalization of the other 
varies while it's spectral shape is unchanged. In this case 
the shape of the frequency resolved spectrum would not 
depend on Fourier frequency and would be identical to 
the spectrum of the variable component. The spectrum 
of the non-variable component could, in principle, be de- 
termined by subtracting the frequency resolved spectrum 
from the average spectrum with appropriate renormal- 
ization. Importantly, in this example the X-ray flux in 
all energy channels will vary coherently and with zero 
time/phase lag between different energies. Presence of sig- 
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Fig. 3. The low (left) and high (right) frequency parts of the power spectrum of 4U1608 averaged over all data used for 
analysis. The power spectrum of the high frequency part, showing two kHz QPO peaks was obtained by "shift-and-add" 
method. 



nificant phase lag and/or Fourier frequency dependence of 
the frequency resolved spectra would indicate that a more 
complex pattern of spectral variability is taking place and 
interpretation is then less obvious. With few exceptions 1 , 
phase lag between light curves in different energy bands 
in luminous LMXBs (e.g. Sco X-l, GX5-1, 4U1608-52, 
4U0614+091 etc.) is usually small, A<j) ;$ few x 10~ 2 
coherence is consistent w i th un ity, (e.g. IVaughan et all 
I1994L Il99flt iDieters et all 1200(1) and the fractional rms- 
vs-e nergy dependence si milar at different Fourier frequen- 
cies ijvan der KlisLll98l . This suggests, that Fourier fre- 
quency resolved spectral analysis can be applied and its 
interpretation is sufficiently straightforward and model- 
independent. 

The structure of the paper is as follows. We briefly de- 
scribe the data in Sect. [21 In Sect. El we present the results 
of the observations, show that the frequency resolved spec- 
tra do not depend upon the Fourier frequency, and con- 
strain the time lags between different energies. The initial 
observational results are summarized in Subsect. 13.31 In 
Sect. 0] we show that a particularly simple form of the 
spectral variability is required in order to satisfy the ob- 
servational constraints - the flux variations in different 
energy channels must be related by a simple linear trans- 
formation. In Sect. El we compare the expected spectra of 
the disk and boundary layer emission with the frequency 
resolved spectra. We show that the observed aperiodic 

1 ~ 150° phase lag was detected in the normal branch QPO 
of Cyg X-2 with the pivot energy ~ 5 keV, bu t no such lags 
were found in a similar spectral state of Sco X-l jDieters et all 
120001) 



and quasiperiodic variability is primarily caused by vari- 
ations of the luminosity of the boundary layer and its en- 
ergy spectrum can be represented, to certain accuracy, by 
the frequency resolved spectra. In Sect. El we discuss the 
boundary layer emission spectrum, it's dependence on the 
mass accretion rate and implications of these results for 
disk and boundary layer models. The results are summa- 
rized in Sect. |3 

2. The data 

In order to investigate the spectral variability of neutron 
star LMXBs in the soft (high luminosity) state we have 
used the RXTE/PCA l|Bradt et all Il993|) observations 
of one of the bright Z-sources (Hasinger & van der Klis 
1989) GX340+0. Our choice was defined by the require- 
ment that the PCA configuration combined sufficiently 
high energy resolution (large number of the energy chan- 
nels) with good timing resolution and large total exposure 
time. We selected observations from proposal P20053 per- 
formed from Sep. 21 to Nov. 4, 1997 with total exposure 
time of «178 ksec. The PCA configuration provided 31 
energy channels in the 3-20 keV energy band with 2 msec 
time resolution. The detailed timing analysis of these data 
was presented earlier by Jonker et al. (2000). 

The intensity-color diagram of GX340+0 is presented 
in Fig.Q] The Fig.EJshows an example of the power density 
spectrum, corresponding to the beginning of the horizon- 
tal branch on the intensity-color diagram. The prominent 
QPO peak at ~ 25 Hz corresponds to so called Horizontal 
Branch Oscillation (HBO). 
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Fig. 4. Average and frequency resolved spectra of 
GX3404-0 in the Horizontal Branch (Z=0-1, see Fig-fl}. 
The solid lines show the Comptonized emission spectrum 
rescaled to match the frequency resolved spectra and the 
high energy part of the average spectrum. The parameters 
of the Comptonized spectrum are best fit to the frequency 
resolved spectrum of the ~ 25 Hz QPO (Table El- 

The kHz QPOs are known to be rather weak in the 
case of GX340+0 (especially the upper kHz QPO), and 
the available data on the source did not have enough sen- 
sitivity to study them in detail. Therefore, to investigate 
the properties of two kHz QPOs and compare them with 
the lower frequency variability we chose another bright 
LMXB, 4U1608-52. The data of proposal P30062 provide 
energy and timing resolution adequate for this purpose. 
Detailed ana lysis of kHz QPOs in thes e observations was 



presented in lMendez et al.1 l)1999ll2001|) . During these ob 
servations the source was in the high luminosity state and 
had two strong kHz QPO peaks at frequencies ~ 600 — 800 
Hz and ~ 900 — 1200 and somewhat weaker QPO at ~45 
Hz (Fig. 12). Our aim is to compare the spectral variability 
corresponding to these three QPO components. 

The reduction of the PCA data was performed with 
the help of standard tasks from FTOOLS /LHEASOFT, 
version 5.1. The spectral approximations were done using 
XSPEC. In all spectral fits the interstellar absorption was 
fixed at NH = 5 • 10 22 cm" 2 and NH = 1 • 10 22 cm" 2 for 
GX340+0 and 4U1608-52 respectively. 



3. Fourier frequency resolved spectra 

3.1. Low frequency continuum and QPO 

As is well known, the X-ray variability properties and pa- 
rameters of Z-sources depend significantly on the spec- 



2 

2 1-5 

*H 1 



1.4 
2 1.2 

0.8 
1.2 



a 
f — i 



1 

0.8 
1.2 

o 

1 

U 

0.8 

1.4 

o 1.2 

^ 1 
0.8 

0.6 

1.4 

o 1-2 

£ 1 

to 

fn 0.8 
0.6 



1 1 1 1 1 1 1 1 _ 
■V 0.03-0. 16 Hz - 


i i i i i 1 i 


i i i i i i | i 
~ 1 , 0.16-0.40 Hz~ 

+ i 


i 


i i i i i i 




i i 

i+j 


i i i i i i 

0.4-0.9 H 


z — 


L + 

1 


i i i i i i 


H 


| 


i i i i i i 

0.9-7.9 Hz- 

. 11-+-, . ++,+>J^~ I ! 


i- "f 

i 


i i i i i i 






i i i i i i 

h ' 7.9-16 Hz 






i i i i i 1 i 


• 




i i i i i 1 i 






1 32-64 Hz 

i i i i 





10 

Energy, keV 

Fig. 5. GX340+0. Ratios of the frequency resolved spectra 
in various Fourier frequency bands to that of the 16-32 Hz 
band, corresponding to the Horizontal branch QPO. The 
data are from the Horizontal Branch of the color-intensity 
diagram (Z=0-1), i.e. the same as in Fig. 21 Note that the 
vertical scale is different in different panels. 



tral state of the source as given by its position on 
the Z-diagram, the highest level of variability being ob- 
served at the beginn ing of the Horizontal Branch (see e.g. 
Ijonker et alJ l)2000|) '). Therefore throughout most of the 
paper we used the data averaged over Horizontal Branch of 
the color-intensity diagram (shown by a polygon in Fig.rj). 
The results for the normal branch are discussed in Sect.EJ 
The power spectrum of GX340+0 on the horizontal branch 
is shown in Fig. 

The Fourier frequency resolved spectra in several fre- 
quency bands corresponding to the band limited contin- 
uum noise component and the ~ 25 Hz QPO are shown in 
Fig. 0] along with the conventional spectrum of the source 
averaged over the same data. The ratios of the frequency 
resolved spectra to that of the QPO are presented in 
Fig. El The spectra and the ratios demonstrate that their 
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Fig. 6. Phase lags for GX340+0 in the Horizontal Branch 
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The energy dependent phase lags were computed in the 
1-32 Hz frequency range, the frequency dependent lags 
are between 3-6.5 keV and 6.5-13 keV energy bands. The 
phase is normalized to 0-1 interval. 



shape depends on the Fourier frequency at low frequencies 
and becomes independent of th e frequency a t f ^ 0.5 Hz. 
Similar behaviour was found by van der Klisl lll986) in the 
case of another bright LMXB - Z-source GX5-1. In par- 
ticular, it was noted, that Horizontal Branch Oscillations 
and the low frequency noise have the same, hard, frac- 
tional rms spectrum and that at lower frequencies vari- 
abilty becomes softer. 

The upper limits on the possible variations of the spec- 
tral shape with Fourier frequency depend on the photon 
energy and Fourier frequency and vary from ~ 5% to 
- 25% (Fig. EJ. Note that in the 0.5 <, f <, 30 Hz fre- 
quency range these upper limits are determined by the sta- 
tistical uncertainties only. At higher frequencies, / i; 30 
Hz, there might be an indication, although not statisti- 
cally significant enough, of a frequency dependence of the 
spectral shape (lower panel in Fig.^J. Another conclusion 
from the data presented in Fig. 0] important for the fol- 
lowing discussion, is that all frequency resolved spectra 
are significantly harder than the average source spectrum. 

The phase lags as function of energy and Fourier fre- 
quency are shown in Fig. No statistically significant 
phase lags were detected in the Horizontal Branch of 
GX340+0 with upper limit of A<fi ~ 10~ 2 , where phase 
<j> is normalized to the interval 0-1 (as opposed to — 2ir). 
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Fig. 7. Upper panel: Average and frequency resolved spec- 
tra of 4U1608-52. The solid lines show Comptonized disk 
spectra rescaled to match the frequency resolved spectra 
and the high energy part of the average spectrum. The 
parameters of the Comptonization spectrum are the best 
fits to the frequency resolved spectrum of the lower kHz 
QPO (Table |2| . Lower panel: Ratios of the spectra shown 
in the upper panel to the spectrum of lower kHz QPO. 



3.2. kHz QPOs 

As was mentioned above, the strength of the kHz QPOs in 
the power spectrum of GX340+0 is insufficient to study 
their frequency resolved spectra in detail. We therefore 
had to use for this purpose a different source - 4U1608-52 
(Fig. |3J). As the frequency of the kHz QPO varied dur- 
ing the RXTE observations we used "shift-and-add" tech- 
nique l)Mendez et all l200lh in order to improve statistics. 
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During the same observations the source showed a lower 
frequency QPO at ^45 Hz as well. 

The Fourier frequency resolved spectra of the two kHz 
QPOs and the HBO along with their ratios are shown 
in Fig. [7| As in the case of GX340+0, all three spectra 
are consistent with each other within the statistical un- 
certainties. The upper limits on possible variations of the 
spectral shape, however, are significantly less constraining, 
- 30 - 50% at best. 

3.3. Summary of the observational results 

The results of the above analysis can be summarized as 
follows: 

1. The shape of the frequency resolved spectra does not 
depend or depends very weakly on the Fourier fre- 
quency. This includes the continuum band limited 
noise component, the Horizontal Branch QPO and 
the kHz QPOs. The constraints on possible Fourier 
frequency dependent variations of the spectral shape 
are rather stringent for the band limited noise and 
Horizontal Branch oscillations (GX340+0) - ^ 5-25% 
in the 0.5 ^5 / ^ 30 Hz frequency range and signifi- 
cantly less restrictive for kHz QPO - ^ 30 — 50%. 

2. The Fourier frequency resolved spectra at all frequen- 
cies are significantly harder than the average spectrum. 

3. No statistically significant phase lags were detected in 
the case of GX340+0 with upper limits of Ac/) < 10~ 2 
in the 0.1 — 30 Hz frequency range and 3 — 17 keV 
energy range. 

4. At low Fourier frequencies, / 0.5 Hz, frequency de- 
pendent variations of the spectral shape become im- 
portant - the frequency resolved spectra tend to be- 
come softer (GX340+0). 

4. Interpretation of the Fourier frequency resolved 
spectra 

We show here that independence of the frequency resolved 
spectra on the Fourier frequency and the smallness of the 
phase lags observed in majority of bright LMXBs (but 
see footnote 0) require a particularly simple form of the 
spectral variability. 

The constancy of the spectral shape with Fourier fre- 
quency implies that the power spectrum P{E, w) can be 
represented as a product of two functions, one of which 
depends on the energy and the other on the frequency 
only. For convenience we write P{E, ui) in the form: 



P{E,u)=S 2 {E) x/V) 



(1) 



where non-negative functions S(E) and f{u>) can be di- 
rectly determined from the frequency resolved spectra. 
The Fourier image of the light curve F(E, t) is: 



F{E,u) = S{E) x f(u) x e l ^ E ^) 



(2) 



In the general case the complex argument 4>{E, to) can 
depend both on Fourier frequency u> and energy E. If the 



phase lags between different energies are negligibly small, 
4> depends on the Fourier frequency only and the Fourier 
image of F(E,t) is: 



F{E,u) = S{E) x f{u) x e l 



(3) 



The light curve F{E,t) can be computed via inverse 
Fourier transform of F{E,uj): 



F(E,t) = / duF(E,w)e i 



S(E) x / dw/He*' 



(4) 



= S(E) x fit) 



An arbitrary function of energy can obviously be added 
to the above expression. Thus, the light curve at energy 
E should satisfy the equation: 



F{E,t) = S {E) + S{E) x f{t) 



(5) 



i.e. the light curves at different energies are related by a 
linear transformation. 

Eq. J5J significantly restricts the pattern of the spectral 
variability. Suppose that the spectrum of the source at any 
given moment of time can be represented as: 



F(E,t) =S (E)+A(t) xS(E,p(t)) 



(6) 



where A(t) is varying normalization, and p(t) represents 
variations of the spectral parameter, on which the spectral 
flux at a given energy depends non-linearly (e.g. tempera- 
ture, Thompson optical depth etc.). Using Taylor expan- 
sion: 

F(E, t) = S {E) + A{t) x S{E,p{t)) = 
= S {E)+A{t)xS{E,p Q ) + 
dS{E, Po ) 



-A(t) x 



dp 



-A(t) x 



l d 2 S(E, Po ) 
2 dp 2 



x (p(t) - p ) + 



x (p(t) - Po) 



(7) 



In order to satisfy Eq. JSJ, one of the following two con- 
ditions should be fulfilled: (i) p{t) = Po = const, i.e. only 
the normalization varies with time, or (ii) the normaliza- 
tion is constant and variations of the spectral parameter 
are sufficiently small, so that the terms higher than linear 
in the Taylor expansion can be neglected. Simultaneous 
variations of comparable amplitude of the normalization 
and the spectral parameter (or of any two spectral param- 
eters) would be consistent with observations only in the 
case of A(t) x {pit) — Po) — const or in the case of specific 
shape of the power density spectrum of A{t) (e.g. flat). 

From Eq. JSJ it follows that the coherence of the signals 
in any two energy bands is exactly unity, as they are re- 
lated by a linear transformation. This is in agreement with 
the observed behaviour - as illustrated by an example of 
GX340+0, shown in Fig. [SJ the coherence between the 
light curves in 3- 6.5 and 6.5-13 keV energy bands is close 
to un ity (see also lVaughan et al.ll994lll999HDieters et alJ 
2000' for results for other luminous LMXBs). 
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Fig. 8. GX340+0: Coherence between the light curves in 
the 3-6.5 and 6.5-13 keV energy bands as function of fre- 
quency. No correction for the dead time effects has been 
made. 

4.1. Case of black hole binaries 

The spectral variability in the black hole binaries is more 
complicated than described by Eq. J5J because of the 
strong dependence of the shape of their frequency resolved 
spect ra upon the Fourier frequency l|Revnivtsev et all 
Il999lh 2 rather than non-zero time/phase lags observed. 
Indeed, typical values of the phase lags measured for Cyg 
X-l are negligible in this context, A<j> ^ few x 10 -2 (e.g 
iKotov. Churazov fc GilfanovLl200l|) . Therefore 

<j>{E,uj) = 0oM + <M£» 

where 

<pi{E,oj) « 0oM 

and Eq. Q and therefore Eq. © will hold with sufficient 
accuracy if condition of Eq. is fulfilled. 

5. Application to bright LMXBs 

5.1. Frequency resolved spectra and nature of the 
variable component 

Thus, observed properties of the luminous LMXBs 
(Sect. require a particularly simple form of the spec- 
tral variability on time scales from ~ sec to ~ msec, de- 
scribed by Eq. {JJ. The term Sq(E) in Eq. 10) represents 
constant (non-variable) part of the source spectrum and 
f(t) represents either (i) variations of t he normalization 
of th e variable spectral component (e.g. IChurazov et al j , 
12001) or (ii) small varia tions of the spectral parameter 
fe.g.lMendez eTaTHTflfl^. 

If variability of the X-ray flux is caused mainly by vari- 
ations of the normalization, it follows from Eq. JQ), that 
the spectrum of the variable component S(E) is identical 

2 With a possib le exception of the sof t state, e.g. in the soft 
state of Cyg X-l llChurazov et al.ll200ll) found that the shape 
of the power density spectra does not depend on energy and 
that the spectral variability satisfies the Eq. @. 



in shape to the frequency resolved spectrum, i.e. can be 
directly determined from observations. The spectrum of 
the non- variable component can, in principle, be obtained 
subtracting S(E) from the averaged spectrum with appro- 
priate renormalization. 

In the second case (small variations of the spectral pa- 
rameter) , the energy dependence of the frequency resolved 
spectrum S(E) corresponds to the first derivative of the 
spectrum with respect to the parameter p : dS< ^ ) ' p ^ , which 
might differ significantly from the spectrum itself. Indeed, 
considering a Wien spectrum , S(E T) = E 2 e~ E / T , as 
an example (cf. iMendez et alJ l)l997(l h the intensity vari- 
ations are: 



F(E,t) = S(E,T ) + 
= S (E) 



dS(E,T )) 
dT 



x ST(t) 



(8) 



, N E ST(t) 
So(E) x - x — ii 
Jo Jo 



Therefore in Eq. © S(E) = S {E) x E/T , i.e. is harder 
than the average spectrum. It is intuitively obvious, how- 
ever, that in order the linear expansion to be valid, the 
variations amplitude at any energy should be smaller than 
the average flux. Indeed, for the quadratic term in Eq. 
to be neglectable, the following condition should be satis- 
fied: 



E ST(t) 
— x — — 

T T 



< 1 



(9) 



i.e. (cf. Eq. JSJl) the frequency resolved spectrum in any 
frequency range and at any energy should be small in com- 
parison with the average spectrum of the variable compo- 
nent. Equivalcntly, fractional rms of flux variations com- 
puted with respect to the average flux of the variable com- 
ponent (as opposite to the total average flux) should be 
less than 1. 

In the first case identification of the variable compo- 
nent would be easy and unambiguous. However, no dis- 
tinction between the two possibilities outlined above can 
be made based solely on the results of the frequency re- 
solved analysis, and additional considerations should be 
taken into account. 

5.2. Emission from the boundary layer and accretion 
disk 



Base d on theo re tical g rounds llSunvaev fc Shakural 
I l986t iKluzniaki Il98s( hnogamov SunvaevL I199PJ: 
IPopham fc Sunvaev . l200ll) and observational r esults 
fe.g. iMitsuda et all Il984t iMitsuda k TanakaL Il986t 
iGierlinski fc DoneL I 20oJT it is expected that in the case 
of accretion onto a slowly rotating weakly magnetized 
neutron star at sufficiently high mass accretion rates 
there are two major components of the accretion flow: 

1 . the optically thick accretion disk extending close to the 
surface of the neutron star or last marginally stable 
orbit 
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Fig. 9. The average and frequency resolved spectra of GX340+0 (left) and 4U1608-52 (right). The GX340+0 data are 
from horizontal branch of the color-intensity diagram, i.e. the same as used before. The shaded area shows a plausible 
range of the boundary layer spectral shape calculated subtracting the predicted disk spectrum from the total spectrum 
and renormalizing the residual to the total energy flux of the frequency resolved spectrum (see Sect. l5.2l and Tablc^for 
details). The dashed histogram shows the accretion disk spectrum with parameters from Tablel2"lfSect. l5~3")) . the upper 
solid histogram shows the difference between the total and accretion disk spectrum (« boundary layer spectrum). The 
lower solid histogram is the same but scaled to the total energy flux of the frequency resolved spectrum. The lower 
panels show ratio of the frequency resolved spectrum to the lower histogram. 



2. the boundary la yer (e.g.|Pojahain_^_Si inYa£ or 
spreading layer l|Inogamov fc Sunvaevi Il999» near the 
surface of the neutron star. In this layer the accreting 
matter decelerates to the spin frequency of the star 
and spreads over its surface. 

Simple theoretical arguments, taking into account energy 
released in the accretion disk and in the boundary layer, 
and difference in their emitting areas suggest that the 
spectrum of the boundary layer emission should be no- 
ticeably hard er than that of the optically thick accre- 
tiondisk fe.g.lMitsuda et all Il983 iGrebenev fc Sunvaevi. 
2002). Considering the restrictions on the character of 
spectral variability imposed by Eq. JSJ, it is natural to 
assume, that the observed variations of the X-ray flux 
should originate in one of these components. Furthermore, 
X-ray variability should be mainly caused by either varia- 
tions of the total luminosity with nearly constant spectral 
shape or by small variations of the spectral shape. In order 
to distinguish between these possibilities and to identify 
the nature of the varying component we consider below 
theoretical expectations for the disk and boundary layer 



spectra and compare them with the observed frequency 
resolved spectra. 

Due to complexity of the boundary/spreading layer 
problem, no models capable to directly predict its emis- 
sion spectrum exist yet. An attempt of solve the ra- 
diation transfer problem using density and tempera- 
ture profil es obtained in the hydrod ynamical simula- 
tions (e.g. IPopham fe Sunvaevi feoOlj) failed to repro- 
duce, even qualitatively, the observed X-ray spectra and 
demonstrated importance of the self-consistent treatment 
llCrebenev fc Sunvaevi Eool. 

Significantly better progre ss has been achieved with 
the spectra of accretion disks dShakura fc Sunvaevi Il973t 
IShimura fc Takaharal Il99i iRoss fc Fabianll 199fijT It has 
been shown that at sufficiently high mass accretion 
rates the effects of Compton scattering can be approx- 
imately accounted for by introducing a dilution factor 
to describe deviation of the spectr a emitted at each ra- 
dius of the disk from black body (Sh imura fc Takaharal 
Il995l IRoss fc Fabianl \ 199fih . Additional modifications of 
the disk spectrum arise due to gravitational redshift 
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Table 1. Range of parameters used to estimate plausible 
range of the accretion disk and boundary layer spectra 
(Fig.© 



Parameter 


GX340+0 4U1608-52 


D, kpc 


8.5-10.5 3.5-4.5 


-Fb.i-30 kcv 1 


2.7 ■ 10~ 8 8.8 ■ 10~ 9 


cos(i) 


0.3-0.7 


Tool/Toff 


1.6-2.0 


v NS ,Hz 


0-700 


Vtot i 


0.213-0.116 


/disk 


0.26-0.57 


f 2 

JBL 


0.74-0.43 



1 - absorption corrected, erg/s/cm 2 

2 - derived values of the total accretion efficiency and 
fraction of the energy released in the disk and boundary layer 
respectively, corresponding to the considered range of the 
neutr on star spin frequency, for a 1.4Mq neutron star, EOS 
FPS (|Sibgatullin fc Sunvaevil2000h . 



and D o ppler effects llEbisawa, et all 199lllRoss fe Fabianl 
I 1996t Bhattacharvva et all l200llh Relatively simple 
models of multicolor disk type which account for the 
effects of Compton scattering with a simple color-to- 
effective temperature ratio turned out to be successful 
in describing the accretion disk sp ectra observed in the 
high state of black hole systems fe.g. Ebisawa et all ll991: 
lOierliuski et all I 1997h. 

Therefore we use the predicted disk spectrum as i\ 
starting point and subtract it from the total spectrum 
in order to calculate expected spectrum of the boundary 
layer. To estimate the plausible range of the boundary 
layer spectra we investigate the parameter space of the 
accretion disk model. This approach enables one to pre- 
dict the disk and the boundary layer spectra based on 
the observed X-ray flux and spectrum and very generic 
system parameters, such as neutron star spin frequency 
or the source distance. For this analysis we adopt the 
genera l relativistic accretion disk model bv lEbisawa et all 
l|l99l|) (the "grad model" in XSPEC). The parameters of 
the model are: the source distance D, mass of the cen- 
tral object Mns, disk inclination angle i, the mass accre- 
tion rate M and the color-to-cffectivc temperature ratio 
/ = T co i/T e ff. The source distance was varied in the range 
8.5-10.5 kpc (GX340+0) and 3.5-4.5 kpc (4U1608-52), 
the disk inclination - in the range of cosi = 0.3 — 0.7, 
color-to-effective temperatur e ratio - in the range 1.6-2.0 
l|Shimura fc Takaharai [1995) . 

In order to estimate the plausible range of mass 
accretion rates we used the total accretion efficiency 
rytot and relative energy release fractions of the disk, 
friisk, and the neutron star surface, Jbl, calculated by 
Sibgatul iinfc Sunvaevl ll2000lh These account exactly for 
the space-time metric around a neutron star with a given 
spin frequency and equation of state. In particular we used 
their results for a neutron star with gravitational mass of 
1.4 Mq described by FPS equation of state (their eqs. 



(1) and (2)). In computing the relation between the ob- 
served luminosity and the mass accretion rate we ignored 
the light bending and abberation effects, the real geome- 
try of the spreading layer and shadow cast by the neutron 
star on the disk and boundary layer: 

L ohs = (d(i)f disk + d(| - i)f B h) x 7/ to t Mc 2 (10) 

where emission diagram is described by Chandrasekhar- 
Sobolev law d(i) ss 3/7(1 + 2. 06 cosi) cosi. The luminosity 
iobs was computed from the absorption corrected fluxes 
extrapolated to 0.1-30 keV energy range. The total cor- 
rection factor from observed 3-20 keV flux was ~ 2 in 
both cases. Therefore, unless a new spectral component 
appears outside the 3-20 keV PCA energy range, the un- 
certainty of the flux correction should not exceed few tens 
of per cent. The accretion efficiency r^ot and the disk and 
boundary layer fractions /disk and /bl depend, obviously, 
on the neutron star spin, which was varied in the range 0- 
700 Hz assuming that the neutron star and disk co-rotate. 

The investigated range of the disk parameters is sum- 
marized in Table ^ Each predicted disk spectrum was 
subtracted from the total spectrum and the residual 
(^boundary layer spectrum) was normalized to the ob- 
served 3-20 keV energy flux of the frequency resolved spec- 
trum. The resulting range of the boundary layer spectra 
for two sources is shown in Fig. as the shaded area along 
with the average and frequency resolved spectra of the 25 
Hz and kHz QPO and the best fit model described in 
Subsect. 15.31 The lower panels in Fig. [5| show ratio of the 
frequency resolved spectrum to predicted boundary layer 
spectrum. Different models of the disk emission, e.g. ignor- 
ing the relativistic effe cts (the "dis kpn" model in XSPEC , 
iGierlinski etldll 1997|) or model of lRoss fc Fabianl l) 1996|) 
which more accurately accounts for the Compton scat- 
tering effects, but assuming Newtonian dynamics, change 
the details but do not alter the general conclusion. The 
close similarity of the expected boundary layer and the 
observed frequency resolved spectra strongly suggests that 
X-ray variability is caused by variations of the luminosity 
of the boundary layer component. It also implies that the 
spectrum of the boundary layer emission does not change 
significantly in the course of these variations (cf. example 
given by Eq. ©). 

On the other hand it is obvious from Fig. that the 
disk spectrum is significantly softer and can not reproduce 
the observed frequency resolved spectrum. Qualitatively, 
this fact does not rule out the possibility that X-ray vari- 
ability is due to variations of the disk emission spectrum. 
Indeed, variations of e.g. disk temperature can result in 
the frequency resolved spectra significantly different (in 
particular, harder, cf. Eq. ©) than the average spectrum. 
However, as discussed in Sect. 15.11 in order to maintain 
constancy of the shape of the frequency resolved spectra 
and zero time lags between different energy channels it 
is required, that frequency resolved spectra in any fre- 
quency range are <C than the average spectrum of the 
variable component. This condition is obviously violated 
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Table 2. Parameters of the spectral modeling 



4U1608-52 
kHz QPO 



10 - 



10 

Energy, keV 

Fig. 10. The energy dependence of the fractional rms of 
the kHz QPO in 4U1608-52. Flattening above ~ 10 keV 
indicates, that the energy spectrum of the oscillations at 
these energies has the same shape as the total spectrum. 



at E <; 10 — 15 keV (Fig.[§|- cf. the disk and the frequency 
resolved spectra). 

Considering the quick decline of the disk emission 
above ~ 10 keV, another test of the suggested interpre- 
tation might be the behavior of the frequency resolved 
spectrum (or fractional rms) at high energy, <; 10—15 
keV. Indeed, as at E J; 10 — 15 keV the total spectrum is 
dominated by the boundary layer emission, the fractional 
rms in the above scenario is expected to become inde- 
pendent of the photon energy. Due to insufficient energy 
coverage this can not be checked in the case of GX340+0 
where the disk temperature is rather high. Flattening of 
the fractional rms is clearly seen in 4U1608-52 (lower-most 
panel in Fig. Hand Fig. llO|) . having ~ 10 times lower mass 
accretion rate and correspondingly lower disk temperature 
an d softer disk spectrum . Similar behavior was observed 
bv lMitsuda et all l)l984|) for Sco X-l using TENMA data. 

Thus, we conclude that the bulk of the X-ray variabil- 
ity observed in GX340+0 and 4U1608-52 and, presum- 
ably, in other luminous LMXBs, on ~ second - millisecond 
time scales is due to variations of the boundary layer lumi- 
nosity. The shape of the boundary layer emission spectrum 
remains nearly constant in the course of these variations. 
Therefore the energy spectrum of the variable component 
- the frequency resolved spectrum should be representa- 
tive, to some accuracy, of the spectrum of the boundary 
layer emission. This can be used to separate the bound- 
ary layer and the accretion disk contribution to the total 
spectrum and permits to check quantitatively the predic- 
tions the accretion disk and boundary layer models. It also 



Parameter 



GX340+0 



4U1608-52 



D, kpc 
NH, cm" 2 

F 3 _2o kev 1 erg/s/cm 2 
Fo.1-30 kev 2 erg/s/cm 2 
£o.i-30k c v 2 erg/s 
M 3 , g/s 



8.5 
5 ■ 10 22 
1.4 • 10" 8 
2.7- 10~ 8 
2.3 • 10 38 
3.1 • 10 18 



4.0 
1 ■ 10 22 
4.2 ■ 10~ 9 
8.8 ■ 10" 9 
1.7 ■ 10 37 
2.0 ■ 10 17 



QPO frequency resolved spectra (^boundary layer) 
power law with exponential cutoff (phabsxcutoffpl) 



-0.55 ±0.16 



-1.28 ±0.13 



E f , keV 


3.3 ±0.2 


2.4 ±0.1 


X 2 /dof 


13.9/16 


4.8/9 


Comptonization model (phabsxcomptt) 


kT hb , keV 


1.3 ± 0.2 


1.4 ±0.4 


kT c , keV 
r 

X 2 /dof 


o i +0.9 
o,x -0.S 

6.o±|:f 

11.3/15 


2.6t °°3 
4.0/8 


average spectra 4 


phabs x (grad+pexrav±gaussian) 


inclination 5 60° 


70° 


Tcol/Teff 


1.7 


1.8 


M, 10 18 g/ 


s 3.0 ±0.04 


0.34 ±0.01 


M/M Edd 6 


w 0.9 


» 0.1 


Q,/2n 


0.27 ±0.07 


£ 0.1 


EW 7 , eV 


48 ± 10 


147 ± 17 


LBL/L to t, 3-20 keV 47% 


57% 



1 - observed; 2 - absorption corrected; 3 - calculated 
from the total unabsorbed luminosity using the accretion 
efficiency for the neutron star spin frequency of z/ns = 500 Hz 
iSibeatullin fc Sunvaevl l200l see Sect. EH Eq. (JTOJl); 4 - the 
details of spectral modeling are given in Sect. 15.31 5 - fixed at 
fudicial value; 6 - assuming MEdd = 2- 10 38 /c 2 ?7di s k ~ 3.5- 10 18 
g/s, where 77disk = 0.066 - accretion disk efficiency for 
vns = 500 Hz and Mns = IAMq; 7 - line energy and width 
were fixed at Fii ne = 6.7 keV and an ne = 0.5 keV 



opens the possibility to measure relative contributions of 
these two components of the accretion flow to the total 
observed X-ray emission. 

5.3. Spectral modeling 

As small variations of the spectral shape of the boundary 
layer emission in the course of flux variations can not be 
excluded and are likely to take place, the frequency re- 
solved spectra represent the boundary layer emission with 
a certain accuracy only. This should be kept in mind while 
interpreting the spectral fitting results described in this 
subsection. However, the reasonable quality of the spectral 
fits and values of the spectral parameters for two sources 
with significantly different mass accretion rate provide ad- 
ditional indirect support to the main conclusion of this 
section. 

Below we approximate the frequency resolved spec- 
trum (» spectrum of the boundary layer emission) with a 
plausible model and then fit the total spectrum adding the 
accretion disk component. This procedure is equivalent to 



Gilfanov et al.: Boundary layer, accretion disk and variability in LMXBs 



11 



subtracting the renormalized frequency resolved spectrum 
from the total spectrum, with the renormalization coeffi- 
cient determined by requirement to minimize \ 2 ■ 

The simplest models of the blackbody or Wien spec- 
trum with kT ~ 2 — 2.2 keV, although they repro- 
duce approximately the shape of the frequency resolved 
spectra, show statistically significant deviations from 
the data (Fig. ITT|) and give unacceptable value of \ 2 ■ 
Significantly better the spectra can be described by 
the C o mptonized emission mo d els (|Sunvaev & TitarchukL 
Il980t IPoutanen fc SvenssonL Il996|) . Although simple 



Comptonization models are not expected to be directly 
applicable to the emission spectra emerging from the 
boundary layer, we give in Table HI th e best fit para me- 
ters for comptt model in XSPEC HjTitarchukl. Il994j) for 
easy comparison with other results. A convenient parame- 
terization of the frequency resolved spectra is provided by 
a power law spectrum with exponential cut-off, which we 
will use for the following analysis of the average spectra. 

As the second step, we approximate the average spec- 
tra by a two component model including the disk, bound- 
ary layer emission and the reflected component appear- 
ing due to reflection of the latter from the accretion disk: 
phabsx (grad+pexrav+gaussian). The shape (but not the 
normalization) of the boundary layer emission was fixed 
at the values determined from approximation of the fre- 
quency resolved spe ctrum. The reflection was mod eled us- 
ing pexrav model l|Maedziarz fc Zdziarski Il995jl and a 
broad gaussian line at 6.4-6.7 keV. All parameters of the 
grad model were fixed except the mass accretion rate. The 
mass of the neutron star was fixed at IAMq, The disk in- 
clination angle and the color-to-effective temperature ra- 
tio were fixed at somewhat arbitrary, although reasonable, 
values, which were chosen to approximately minimize the 
residuals. Note, that presence of the line at E ~ 6 — 7 
keV in the average spectrum is statistically significatly 
required by the data, especially in the case of 4U1608- 
52. The parameters of the spectral fits are summarized 
in the Table The best fit models and different compo- 
nents are depicted Fig. [5] and ^2 Although the \ 2 values 
of the fits to the average spectra are in the range 3.0-4.5 
per d.o.f. and are formally unacceptable, the models de- 
scribe the observed total spectra (having very statistical 
significance, signal/noise~ 10 2 — 10 3 per energy channel) 
reasonably well, with relative accuracy of 1 — 2% in the 
3-20 keV band. 

6. Discussion 

The best fit values of the mass accretion rates obtained 
from the spectral fits agree, within a factor of ^ 1.7 with 
the values predicted from the observed energy flux and the 
accretion efficiency expected for a neutron star with spin 
frequency of ~ 500 Hz (Eq. JTUJ) . This fact is especially 
encouraging, as the two sources have accretion rates dif- 
ferent by a factor of ~ 10 (TableJ^J. For such difference in 
the accretion rate the expected disk temperatures should 
differ by a factor for ~ 1.7 — 1.8, resulting in the different 
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Fig. 11. The absorption corrected average (upper panel), 
frequency resolved (~ boundary layer, middle) and ac- 
cretion disk (lower) spectra of GX340+0 and 4U1608-52. 
The solid, dashed and dotted lines in the upper panel show 
the best fit total, boundary layer and accretion disk spec- 
tra. The thin solid lines in the middle panel are best fit 
Comptonized spectra, the thick grey line is Wien spec- 
trum with kt = 2.1 keV. The frequency resolved spectra 
are normalized to the same energy flux. The solid lines 
in the lower panel show best fit spectra of the accretion 
disk, the dashed line is 4U1608-52 spectrum multiplied 
by a factor of 12. See Table [21 for values of the spectral 
parameters. 
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shape and total flux of the disk spectra (Fig. ^1 lower 
panel). However, after the contribution of the boundary 
layer is accounted for, the relativistic disk emission model 
is capable of reproducing both spectral shape and normal- 
ization. 

Despite large difference in the mass accretion rate in 
the two sources, the energy spectra of the boundary layer 
emission are very similar to each other. This is in line 
with the finding that variations of the boundary layer 
luminosity in the broad range of time scales from ~sec 
to ^msec are not accompanied by significant variations 
of the spectra l shap e. Si milar behavior was found by 
iMitsuda et all l)l984 and iMitsuda k Tanakal l)l986|) on 
longer time scales of ~ 10 3 sec. 

Due to short light travel time of the accretion disk in 
the vicinity of the neutron star, ~msec, the reflected com- 
ponent, if originating in the inner disk, could contribute to 
the variable emission and cause deviation of the frequency 
resolved spectra from the true boundary layer spectrum. 
This can not be directly verified with the present data - 
upper limit on the equivalent width of the 6.4-6.7 keV line 
is in the case of both sources ss 110 eV (90% confidence). 
However, for the observed shape of the spectrum of the 
boundary layer, contribution of the reflected component 
with Vt/2-K ~ 0.2 — 0.3, if any, would not exceed 10% in the 
~ 10 — 20 keV energy range, i.e. is comparable or smaller 
than other uncertainties involved. 

Further along the Z-track of GX340+0 in the color- 
intensity diagram, on normal and flaring branches, the 
fractional rms of the X-ray variability decreases signifi- 
cantly, by a factor of ~ 5 — 10. However, the statistics is 
sufficient to place meaningful constrains on the behavior 
of the frequency resolved spectra at the first half of the 
normal branch (Fig. QJ. The data indicates that the be- 
havior of the frequency resolved spectra does not change 
its character - at sufficiently high frequency, / > 1 Hz, 
their shape does not depend upon the Fourier frequency 
and is significantly harder than the average spectrum and 
expected spectrum of the accretion disk. Therefore, in the 
same line of arguments as above, it is representative of the 
spectrum of the boundary layer spectrum. Fit to the fre- 
quency resolved spectrum by Comptonization model re- 
quires infinitely large values of the Comptonization pa- 
rameter. Correspondingly, the boundary layer spectrum in 
the normal branch can be well fit by Wien or blackbody 
spectrum (which are close to each other in the 3-20 keV 
range) with the best fit temperature of kT w 2.4 keV. As 
evident from Fig. llll and. especially, from Fig.^|the high 
energy part, E<; 8 keV of the spectrum of 4U1608-52 and 
horizontal branch of GX340+0 also follows Wien spectrum 
with temperature in the range kT ~ 2.1 — 2.3. The com- 
posite fit of the total spectrum with the disk + boundary 
layer spectrum, the same as in Subsect. 15.31 gives a best 
fit value of the mass accretion rate of M rj 4.6 • 10 18 g/s, 
i.e. higher than in the horizontal branch (cf. Table|2J). This 
is consistent with the commonly accepted interpretation 
that the mass accretion rate increases along the Z-track 
on the color-intensity diagram. 
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Fig. 12. The absorption corrected frequency resolved 
spectra of QPO (« boundary layer emission) in GX340 in 
the horizontal branch, Z=0-1 (lower M) and upper half 
of the normal branch, Z=l-1.5 (higher M). See Fig. \\\for 
specification of the regions in the color-intensity diagram. 
The horizontal branch data is same as in Fig.l4ll9landllll 
The solid line shows Wien spectrum with kT = 2.4 keV. 



The frequency resolved spectra (^boundary layer spec- 
tra) on the normal and horizontal branch of the color- 
intensity diagram are plotted in Fig. 1121 along with Wien 
spectrum with kT — 2.4 keV. Combined with the middle 
panel in Fig. 1111 this plot shows trend in the dependence 
of the boundary layer spectrum upon the mass accretion 
rate in the range M ~ (0.1 — 1.0)MEdd- We can tentatively 
conclude that with increase of the mass accretion rate up 
to a value close to critical Eddington rate the boundary 
layer spectrum in the 3-20 keV energy range approaches 
a Wien spectrum. Interestingly, at lower values of M the 
character of the deviations of the boundary layer spec- 
trum from the Wien spectrum is similar to that expected 
in the situation when Compton scatterings are important 
factor of the spectral formation in the media with inhomo - 
geneous temperature distribution (jRoss fc Fa,hianL l~i996y 
In particular they are qualitatively s imilar to the numer- 
ical results of lGrebenev fc Sunvaevl l|2002F) on the forma- 
tion of the spectrum of the boundary layer. 

The relatively weak dependence of the shape of the 
boundary layer spectrum upon the mass accretion rate 
(Fig. ^] and I12f> and the relative constancy of the Wien 
temperature is somewhat surprising. It implies that in 
the considered range of M ^ O.lMEdd the plasma tem- 
perature at Comptonization depth in the boundary layer 
weakly depends upon the mass accretion rate, i.e. increase 
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frequency, Hz 

Fig. 13. The power density spectra of GX340+0 in the 
3-7 keV energy range (horizontal branch, same data as in 
Fig.l4*ll9land llll) . The histogram with the error bars show 
power spectra of total, disk and boundary layer emission, 
obtained using the procedure described in the Sect. El the 
straight solid line shows a power law P(f) = 3.5e — 5 x 
/ . The power density is shown in units of fractional rms 
per Hz and all three spectra are normalized to the total 
count rate in the 3-7 keV band. Note, that the lower- 
most frequency bin is affected by the windowing effects, 
suppressing the power. 



of the M does not change significantly vertical tempera- 
ture structure in the boundary layer. 

The fact that kHz QPO show the same behavior as 
other components of the aperiodic variability indicates, 
that they have the same origin, i.e. are caused by the vari- 
ations of the luminosity of the boundary layer. Although 
the kHz "clock" can be in the disk or due to it's interaction 
with the neutron star, the actual modulat ion of the X- 
ray flu x occurs on the neutron star surface. iMendez et alJ 
suggested similar interpretation of kHz QPO in 
4U0614+09. In particular they showed, that energy spec- 
trum of the kHz QPO can be approximately described by 
a blackbody spectrum with kT ~ 1.5 — 1.6 keV. Note, how- 
ever, that they found different energy depedence of contin- 
uum aperiodic variability at lower frequencies. That can 
possibly be explained by the fact that the source was in 
significantly lower luminosity state, M ~ 10 _2 MEdd and 
it's energy spectrum had a distinct power law component 
which could dominate variability at lower frequencies. 

The disk emission is significantly less variable and does 
not contribute significantly to the variability of the X-ray 
flux at / <; 0. 5 — 1 Hz. This conclu sion is in agreement 
with results of lChurazov et al.1 ( 200 lh for the high state of 



Cyg X-l, indicating that stability of the X-ray emission 
might be a common property of the optically thick accre- 
tion disk, independently on the nature of the compact ob- 
ject. The origin of the variable component, however, is dif- 
ferent in the case of Cyg X-l (and presumably in the soft 
state of other black hole binaries). Indeed, the spectrum 
of the variable component in the soft state of Cyg X-l is 
identical to the time average spectrum of the hard spectral 
component and is adequatly represented by unsaturated 
Comptonization in hot (kT c ~ 50 — 100 keV) and optically 
thin (tt Si 1) coronal flow with possible contribution of 
non-thermal Comptonization. Variable component in lu- 
minous LMXBs considered in this paper, if interpreted in 
the framework of the Comptonization model, requires sat- 
urated Comptonizaion (Comptonization parameter y ~ 1) 
in the relatively low temperature (kT c ~ 2 — 3 keV) plasma 
and is inconsistent both qualitatively and quantitatively 
with the corona models usually applied to black hole bi- 
naries. 

In the bright LMXB systems, the contribution of the 
disk variability becomes noticeable at lower frequencies, 
below ~ 0.5 Hz, where the frequency resolved spectrum 
changes i t's shape and beco mes softer (Fig. 0] and |5j cf. 
results of Ivan der Khsll 19861 for GX5-1). This can be used 
to estimate the contribution of the disk to the observed 
variability of the X-ray flux. The result depends on the 
character of the disk variability. In order to make a crude 
estimate we assume that the disk variations also obey a 
simple linear relation described by Eq. J5J: 



F(e, t) = F disk (E, t) + F BL (E, t) « 
S di sk(£) x /dtak(t) + S BL (E) x / BL (i) 



(11) 



If disk and boundary layer variations were uncorrelated, 
the power density of the total signal would be: 

P(E,u) ex S disk (E) 2 x 1/diskMI 2 + S BL {E) 2 x I/bi»| 2 

where within the accuracy of this consideration one can 
assume that Sdisk(E) and Sbl(E) are the disk and bound- 
ary layer spectra determined in Sect. 15.31 The functions 
|/disk(k-0| 2 and |/bl(w)| 2 after appropriate renormaliza- 
tion represent power density spectra of the disk and 
boundary layer and can be determined from linear fit to 
the square of the frequency resolved spectra in each fre- 
quency interval. The power density spectra thus computed 
are shown in Fig. 1131 along with the total power spectrum 
of GX340+0 in the soft 3-7 keV energy band. The power 
spectrum of the accretion disk flux variations is different 
from that of the boundary layer, does not extend signif- 
icantly to the high frequency domain, and is consistent 
with a power law with slope of —1: 

PusM) « r 1 

The excess power seen at low frequencies, F 0.5 — 1 Hz 
in the soft energy band (cf . Fig. EJ) can be explained as the 
contribution of the disk variations. At higher frequencies 
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the variability is dominated by the boundary layer emis- 
sion, giving primary contribution to quasi-periodic oscil- 
latio ns and the so calle d band limited noise compoinent 
Ce.g. lvan der Klislll986|) . 

7. Summary 

The initial observational results are listed in Subsect. 13.31 
Below we summarize the constrains on the character of 
the spectral variability and implications for the boundary 
layer and accretion disk models. 

7.1. Constrains on the pattern and origin of the 
spectral variability in luminous LMXBs 

1. Using RXTE/PCA observations of two luminous 
low mass X-ray binaries GX340+0 (on the horizon- 
tal/normal branch of the color-intensity diagram) and 
4U1608-52 we show that the shape of the Fourier fre- 
quency resolved spectra on ~ second - millisecond time 
scales does not depend on Fourier frequency (Fig.^0 
and|ZJ. The range of investigated timescales includes 
the band limited continuum noise, the kHz QPO and 
lower frequency QPOs observed at few tens Hz (Fig. [2 
I2J. Combined with the negligibly small phase lags, 
A(p 10~ 2 (Fig. El, this restricts significantly the pos- 
sible pattern of spectral variability of X-ray flux and 
requires linear relation between flux variations at dif- 
ferent energies (Eq. J5J)). 

Considering significant difference in the expected spec- 
tra of the accretion disk and boundary layer the ob- 
served variations should be associated with either one 
of these two major components of the accretion flow. 
The X-ray variability is caused either by variations of 
it's luminosity under constant spectral shape, or by 
small variations of a spectral parameter (e.g tempera- 
ture or optical depth) - Eq. JjJ. 

2. We compared the Fourier frequency resolved spectra 
with the expected spectra of the accretion disk and of 
the boundary layer. The predicted spectra were based 
on the observed energy flux/spectrum and very generic 
system parameters such as the source distance and 
neutron star spin frequency. The frequency resolved 
spectra are well consistent with the range of the bound- 
ary layer spectra expected for plausible range of the 
system parameters (Fig. El Table H). On the other 
hand, they are significantly harder than the expected 
spectrum of the accretion disk. It is unlikely that the 
observed variations are associated with variations of 
the disk luminosity or spectral shape unless the disk 
temperature is ~ 3 — 4 kev, i.e. current accretion disk 
models are inapplicable to the neutron star binaries. 

3. The above suggests that the major part of aperi- 
odic and quasiperiodic variability observed in lumi- 
nous LMXBs above ~ 0.5 Hz is caused by variations 
of the luminosity of the boundary layer. Its spectral 
shape remains nearly constant in the course of the lu- 
minosity variations. This interpretations receives ad- 



ditional support from the constancy of the fractional 
rms with energy at E lb 10 keV, where expected accre- 
tion disk emission vanishes, found in case of 4U1608-52 
(Fig. EH). 

7.2. Implications for the models of the boundary layer 
and disk emission 

The frequency resolved spectrum is representative of the 
energy spectrum of the boundary layer emission. This can 
be used for a more precise decomposition of the spectra of 
luminous LMXBs into accretion disk and boundary layer 
components and for quantitative comparison with predic- 
tions of the theoretical models. In the following we shall 
assume that boundary layer spectrum is identical to the 
frequency resolved spectrum, bearing in mind that this is 
true to certain accuracy. 

1. In the considered range of the mass accretion rate 
M ~ (0.1 — 1)A/Edd, the boundary layer spectrum in 
the 3-20 keV energy range depends weakly on M. Its 
shape is remarkably similar in GX340+0 and 4U1608- 
52 (Fig.llljl. despite the fact that the two sources have 
a factor of ~ 10 difference in the mass accretion rate 
(Table El . 

In the limit of high M, of the order of ~ MEdd (nor- 
mal branch of GX340+0), the boundary layer spec- 
trum in the 3-20 keV energy range can be adequately 
represented by the Wien spectrum with temperature 
kT w 2.4 keV (Fig. El . At lower values of M (4U1608- 
52 and horizontal branch of GX340+0) the spectra are 
better described by Comptonization model with elec- 
tron temperature of ~ 2 — 4 keV and Comptonization 
parameter y ~ 1 (Table 01. Their high energy part, 
E <; 10 keV, is well represented by Wien spectrum 
with temperature of 2.1 — 2.3 keV. 

2. The average spectra can be adequately described by 
the sum of the renormalized frequency resolved spec- 
trum and the accretion disk emission fFig. ITT)). The 
spectrum of the latter is well described by the general 
relativistic accretion disk model. The other parame- 
ters, such as source distance and disk inclination angle 
being fixed at fudicial but plausible values, the best fit 
value of the mass accretion rate coincides, within a fac- 
tor of ^ 1.7 with that inferred from the observed X-ray 
flux and accretion efficiency appropriate for a IAMq 
neutron star with spin frequency of ~ 500 Hz (Table 
I2J. This agreement is especially remarkable, given the 
luminosity and mass accretion rate in the two sources 
differ by the factor of ~ 10. 

3. The accretion disk emission is significantly less variable 
than the boundary layer emission at Fourier frequen- 
cies / <; 0.5 — 1 Hz. The power density spectrum of the 
disk appears to follow a power law Pdisk(/) oc f^ 1 and 
contributes to the overall variability in the soft energy 
band and in the low frequency domain only fFig. I13l) . 

4. The kHz QPOs apear to have the same origin as aperi- 
odic and quasiperiodic variability at lower frequencies. 
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The msec flux modulations originate on the surface of 
the neutron star although the kHz "clock" might re- 
side in the disk or be determined by the disk - neutron 
star interaction. 
5. Finally we point out that in the case of GX340+0 and 
presumably other Z-sources, the above results apply 
to the normal and horizontal branches of the color- 
intensity diagram. The source behaviour on the flaring 
branch, believed to correspond to super-Eddington ac- 
cretion is more complex and is beyond the scope of this 
paper. 
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